Symmetric functions and random partitions 

o 



C/3 



O 



Andrei Okounkov* 



Abstract 



These are notes from the lectures I gave at the NATO ASI "Sym- 



metric Functions 2001" at the Isaac Newton Institute in Cambridge 
' (June 25 - July 6, 2001). Their goal is an informal introduction to 

asymptotic combinatorics related to partitions. 

These notes are based on the three lectures that I gave in in Cambridge 
in July of 2001. Some parts, especially the computations, are given here in 
more detail than in the actual lectures; other parts, such as connections to 
geometry of the moduli spaces of curves, have been omitted. Our goal in 
these lectures was to give an informal introduction to the results contained 
in 123 EH] (there are also many closely related papers such as [HJ EH U21 
■ EH1 EH EH EH] ) , as well as to the general subject of asymptotic combinatorics 

related to partitions. No attempt, however, will be made here to survey this 
already quite large and very rapidly growing subject. Many further references 
to the literature will be pointed out in text. 

I want to take this opportunity to thank many people who made this 
NATO ASI in Cambridge such a wonderful and memorable meeting: the or- 
ganizers, the sponsors, the lecturers, and, of course, the participants. There 
is also a very long list of people who, in addition to their vast influence on 
the field itself, influenced my understanding of it through numerous fruitful 
discussion, in particular: P. Deift, P. Diaconis, S. Fomin, A. Its, K. Johans- 
son, R. Kenyon, R. Stanley, C. Tracy, A. Vershik, H. Widom, and many 
others, including, certainly, my collaborators A. Borodin, G. Olshanski, and 
N. Reshetikhin. 



*Dcpartmcnt of Mathematics, University of California at Berkeley, Evans Hall #3840, 
Berkeley, CA 94720-3840. E-mail: okounkov@math.berkeley.edu. Partial financial support 
by NSF grant DMS-0096246, Sloan foundation, and the Packard foundation is gratefully 
acknowledged. 



1 



It is very tragic that one of the pioneers of the field and, also, one of 
the organizers of this meeting, Sergei Kerov, did not live to participate in it. 
Most of the material in these lectures is directly related to his groundbreaking 
and lasting work and I want to dedicate these lecture to his memory. 

1 Portrait of a large random partition 
1.1 

A partition is a nonincreasing sequence 

A = (Ai > A 2 > A 3 > ••■ > 0) 

of nonnegative integers such that Aj = for all sufficiently large i. Similarly, 
a plane partition is a matrix of nonnegative integers with nonincreasing rows 
and columns and finitely many nonzero entries, for example 

/ 5 3 2 1 \ 
4 3 11 

(1) 

3 2 1 
\2 1 / 

A standard geometric image associated to a partition is its diagram. There 
are several competing traditions of drawing the diagrams of partitions; we 
follow the one illustrated in Figure|H]which shows the diagram of the partition 
(8,5,4,2,2,1). 

Similarly, the diagram of a plane partition is a 3-dimensional object which 
looks like a stack of cubes pushed into a corner. The 3D diagram correspond- 
ing to the partition is shown in Figure ^ 

1.2 

Naturally, partitions are inseparable from symmetric functions. Similarly, 
application of symmetric functions to enumeration of plane partitions is a 
combinatorial classics, an illustrated account of which can be found, for ex- 
ample, in 

Our goal in these lectures is to explain some recent progress in applying 
symmetric function theory to get rather refined information about random 
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Figure 1: A 3-dimensional diagram 7r 



partitions, both ordinary and 3-dimensional. Random partitions, in addition 
to being interesting simply for their own sake, arise in a large variety of other 
contexts. 

One of the most important situations in which random partitions arise 
is the study of increasing and decreasing subsequences in a random permu- 
tation, see for example PQEj- Namely, the celebrated Robinson-Schensted 
correspondence associates to any permutation g 6 S(N) a partition Xns(g) 
of N, which by Greene-Fomin extension (see e.g. [7j) of Schensted's theorem 
[32] encodes the increasing and decreasing subsequences of g. Concretely, the 
number 

Ars(sOi H h X RS (g) k 

of squares in the first k rows of Xns(g) equals the maximal cardinality of a 
union of k increasing subsequences of g. The same is true about columns of 
Xrs and decreasing subsequences of g. 

Increasing subsequences of random permutations play a central role, for 
example, in the analysis of many random growth processes. Therefore, one 
wants to understand the push-forward of the uniform probability measure on 
S(N) under the map g i— > \ns(g)- This push- forward is well known to be 
the Plancherel measure on partitions of N given by 

Ppl(A)= ^, |A| = JV, 
where dim A is the dimension of the corresponding irreducible representation 
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of the symmetric group. 

The number dim A has a very transparent combinatorial meaning: it is 
the number of chains of partitions (also known as standard tableaux) of the 
form 

A = A(°)\A«\A( 2 >\---\A( iV ) = 0, (2) 

where /i \ v means that v is obtained from /i by removing a square. There 
are several useful formulas for the number dim A, in particular, the hook 
formula. The textbooks p?Tl I34*j are a wonderful place to learn about these 
things. 

1.3 

The following law of large numbers for the Plancherel measure was discovered 
by Logan and Shepp [23] and, independently, Vershik and Kerov [3H1 155] - 
Consider a diagram of a partition of N as in Figure El and scale it in both 
directions by factor of (note that the rescaled diagram has unit area). 
The rescaled diagram for a Plancherel-typical partition of a large number 
iV will then look something like the diagram in Figure El Now look at the 
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Figure 2: The limit shape of a Plancherel-typical diagram. 

zigzag path in Figure El formed by the boundary of the diagram. This zigzag 
is a graph of a piecewise linear function and, in this way, the Plancherel 
measure becomes a measure on continuous (and even Lipschitz with constant 
1) functions. The discovery of Logan-Shepp-Vershik-Kerov is that, as iV — > 
oo, this measure converges (see [2H1 EE1 Ell for precise statements) to the 
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delta measure supported on the following function 



~ (warcsin(w/2) + a/4 - u 2 ) , \u\ < 2 , 
\u\ , \u\ > 2 . 



This formula looks, at first, rather formidable, but its derivative is much 
simpler, namely 

Q'(u) = — arcsin — , (3) 

7T 2 

naturally extended by ±1 for u [—2, 2]. 

In fact, the derivative Q'(u) is a very natural object because it measures 
how often the zigzag path in Figure 121 goes up (respectively, down) in the 
neighborhood of a given point on the boundary of the limit shape. In other 
words, if we think of the zigzag path as a random sequence of "ups" and 
"downs", Q'(u) is related to the probability to observe an "up" in a given 
position. Traditionally, probabilities of this kind are called correlation func- 
tions. More precisely, the the probability to observe "up" in a given position 
is the 1-point correlation function of our random sequence of "ups" and 
"downs" . 



1.4 

It is now natural to ask about higher correlations, for example, what is the 
probability to observe several consecutive "ups" , or, more generally, any given 
pattern of "ups" and "downs" ? 

Such questions can be conveniently formalized using the following set 



6(A) = {A,-i + |} C Z + f . 



(4) 



This is the set of "downs" of the boundary of A as illustrated in Figure El 

We remark that the map A i— > 6(A) is a bijection of the set of partitions 
and the set of subsets S C Z + ^ such that 



< oo . 



In other words, the set 6(A) has equally many positive elements and negative 
"holes" (this number is the length of the diagonal in the diagram A). 
We now take a finite set 



X = X(N) = {x u ...,x n }cZ 



i 

2 ' 
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Figure 3: Geometric meaning of 6(A) = {•} 



which is allowed to depend on N in such a way that 

-p= -> Ui e [-2, 2] , N -> oo , 




and all differences X{ — Xj have a limit, finite or infinite. The set X is the 
set where we want to observe the "downs" . The parameters Ui describes the 
global position of our observations on the limit shape Q while the differences 
Xi — Xj record their relative local separation. 

It is clear the following probabilities (a.k.a. correlation functions) 



to observe "downs" in positions X = {xi, . . . ,x n } include all possible infor- 
mation about the local shape of the random partition A. As it turns out, all 
these probabilities have a very pretty N — ► oo limit, namely, we have the 
following: 

Theorem 1 ([5]). For X as above we have 



Ppi (Xc6(A)) 



I 

AH 



E ( dimA ) 2 



xce(\) 



Ppi (X C 6(A)) -> det K sin (xi - x 3 ; fa) 



J i,j=l...n 



(5) 



as N 



oo, where 



K sin (x; 0) 



sin (f)x 



(6) 



TlX 



is the discrete sine kernel with parameters 



4>i = arccos(wj/2) 



i = 



1,... 



n . 



(7) 
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In the formula J6J), it is understood that 



K sin (oo; d>) = 



which implies that the determinant in (j5J) factors into blocks corresponding 
to distinct values of U{. In other words, it means that the occurrences of 
"downs" in different parts of the limit shape become independent in the 
N —>■ oo limit, which, intuitively, is what one expects. 
It is also understood in © that 

K sin (O;0) = ^. 

7T 

In particular, the 1-point correlation function, that is, the density of "downs" 
is 

™ , r , A ,„„ arccosw/2 x 

Ppi ({x} c 6(A)) '- , 7= ¥ u . 



This translates into (jSJ) and also explains the role of the parameter (f> in the 
sine kernel: it is the density parameter. 

Finally, it is natural to extend the function (JJJ) by or it for u ^ [—2, 2]. 
Theorem Q remains valid and says simply that for u > 2 (resp. u < —2) the 
zigzag in Figure 121 goes only up (resp. down). 



1.5 



One encounters a formula which looks just like the formula (0) in the random 
matrix theory, see e.g. ^J|2H1- Concretely, consider a Gaussian probability 
measure Ph on the real vector space of Hermitian N x N matrices H with 
density e~ trH . The spectrum of cr(H) of a random matrix if is a random 
iV-element subset of R. Consider the probability density pu(X) to have 
X C cr(H) for some fixed X = {x\, . . . ,x n }, that is, the probability density 
to observe eigenvalues of H in given positions. Formally, this probability 
density is defined as follows 



p H {xi, ■ ■ ■ , x n ) dxi ■ ■ ■ dx n = P H 



a(H) n [J[xi,Xi + dxi 



n 



In other words, ph(xi, ■ ■ ■ > x n ) dx\ ■ ■ ■ dx n is actually a measure whose values 
on a given set A C K n is related to the probability of finding n eigenvalues 
of H is A. 
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It is important to remember that since pu(X) is a density, it transforms 
nontrivially under a change of variables and, in particular, it scales nontriv- 
ially under the scaling of variables. 

The 1-point correlation function pn(xi) describes the density of eigen- 
values of H. The analog of the limit shape Q in this case is the Wigner 
semicircle law 



V4 



u 2 x 



Ph(x)^— , -= -> u e [-2, 2] , N^oo. 



N 27r JN 



In other words, scaled by v^, the density of eigenvalues converges to the 
semicircle with diameter [—2, 2}. The normalization h= is related to the fact 
that there are N eigenvalues of H and, because they all are concentrated on 
the segment [—2vN,2y/W\, the density of eigenvalues is expected to be of 
order y/~N . 

By the same token, the typical spacing between the eigenvalues is of order 
1/yN. It is then meaningful to ask about the asymptotics of p(X) as 



x i 



u 



e [-2, 2] , V N(xi - xj) -> yt - yj , 



'N 

for some fixed yi G R. One finds (see j!4[ 123] ) that in this limit 



P{X) 



det 



KUyi-yj;r(u))], r(u) = y/l - {u/2f . (8) 



This looks almost identical to (0), with one minor and one not so minor 
difference. The minor difference is a different dependence of the density 
parameter in the sine kernel on the global position u. A more important 
difference is that in © the variables Xi — Xj take only discrete values, whereas 
in (jSJ) we deal with continuous variables. 

In the continuous case, we can get rid of the density parameter in the sine 
kernel by a simple rescaling. After this transformation, the global question of 
eigenvalue density and the local question of spacings between the eigenvalues 
become completely decoupled. Much work has been done in the random 
matrix theory to show that the eigenvalue spacing have the same sine kernel 
form for many much more general measures on the space of matrices, in other 
words, to show that the occurrence of the sine kernel is universal, see e.g. 

oim. 

Comparing, (jSJ) with (JSJ) reveals an striking "universality of formulas", 
namely that the formulas are the same even when the meaning is different. 
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This parallel, however, is more than purely formal. For example, the discrete 
sine kernel process shares many features of its continuous counterpart, such 
as the repulsion of eigenvalues which we will discuss momentarily. Also, 
continuous distributions, such as the distributions of eigenvalues of a N x N 
Gaussian Hermitian matrix, can be obtained as a suitable limit of certain 
natural measures (of the kind considered below) on partitions with N parts. 



1.6 



It can be inferred from the formula (jHJ) that the eigenvalues of a random 
matrix have a strong tendency to be equally spaced, much more so than 
simply a random collection of points on a line. For example, the 2-point 
correlation function (in the unit density case) is equal to 



det 



1 



sin n(y—x) 
Tr(y-x) 



sin ir(x-y) 
n(x-y) 

1 



( si 



sin nix 



y) 



V n(x-y) 

which is very small (has a double zero) near x — y = and has maxima for 



x-y = 1,2,3, 



It was noted by Kerov that a there exists a similar repulsion of "ups" and 
"downs" in Figure El Namely, the "ups" and "downs" tend to be much more 
uniformly spaced than outcomes of independent trials, such as, for example, 
the occurrences of "heads" in independent coin tosses. Let us illustrate this 
by a numerical example. Suppose that we are in the very middle of the limit 
shape (that is, u — 0), where the density of both "ups" and "downs" is equal 
to 



arccos(w/2) 



7T 



u=0 



1 

2 



In this case the discrete sine kernel takes an especially simple form 

sin(7ra;/2) 



[x 



' 2- 



7TX 



1 

2 ' 

(-!)" 

7TX 



0, 



X 
X 
X 



0, 

2k + 1, 
±2, ±4, 



Consider the following probabilities 



a n = P sin (up, down, up, down,. . . , up, down) , 
b n = P sin (up, up,. . . , up, down, down,. . . ,down) 
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where each string has length 2n. These discrete sine kernel probabilities are 
given by certain (2n) x (2n) determinants. If the "ups" and "downs" were 
like "heads" and "tails" of independent coin tosses, both a n and b n would 
equal 2~ 2n . But as the plot of log 2 a n and log 2 b n in Figure |U shows 

a n > 2- 2n > b n , 

which means that the discrete sine kernel process very strongly prefers even 
spacing of its "ups" and "down" . 




2 3 4 5 6 



Figure 4: The plot (from top to bottom) of — log 2 b n , 2n, and — log 2 a. 



1.7 

The connection between the Plancherel measure and random matrices be- 
comes exact near the edge of the Logan- Shepp-Vershik-Kerov limit shape 
and the Wigner semicircle, respectively In other words, the maximal parts 
Ai, A2, • • • of a Plancherel random partition, suitably scaled, behave in N — > 
00 limit exactly like the maximal eigenvalues of an Hermitian random ma- 
trix. For Ai it was proven by Baik, Deift, and Johansson in j2], who also 
conjectured that the same holds for any Aj, i = 2,3,4, .. . (for i = 2 they 
verified this in |3J. 

By Schensted's theorem, mentioned in Section 11.21 this means that the 
longest increasing subsequence of a random permutation is, asymptotically, 
distributed like the largest eigenvalue of a Hermitian random matrix. This 
distribution is known as the Tracy- Widom distribution (33] and is given ex- 
plicitly in term of a certain solution of the Painleve II equation. 
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The random point process describing the maximal eigenvalues of a ran- 
dom Hermitian matrix H is known as the Airy ensemble because the correla- 
tion function of this process are described by the same determinantal formula 
as in (jSJ) with the sine kernel replaced by the Airy kernel 

w«,») = Ai(j)A! ' (,) ~ A! ' (j)Ai(v) - o» 

x y 

Here Ai(x) is the Airy function 

Ai(ar) = — / e is3/3+ixs ds = - cos - + sx ) ds . (10) 

27r J~oo 7T Jo V 3 / 

The joint distributions of the largest, 2nd largest, and so on points of the 
Airy ensemble satisfy quite complicated nonlinear PDEs, generalizing the 
Painleve II equation for the Tracy- Widom distribution. The description of 
the Airy ensemble by its correlation functions is much simpler and this is the 
description that we will use in what follows. 



1.8 

The first proof of the full Baik-Deift- Johansson conjecture was given in [27J 
by establishing a direct asymptotics correspondence between permutations 
and matrices via a certain correspondence between counting maps on surfaces 
and counting branched covering of the sphere. A different proof was soon 
after found, independently, in [5] and [T7j. It is this second proof that will 
be explained in these lectures. 

It is clear from Theorem that all correlation become trivial at u — 2, 
which is because near u = 2 the "downs" become so rare that they occur 
infinitely far away from each other. In other words, for a Plancherel typical 
partition A^ — \j — ► oo as N — * oo for any fixed i ^ j. In fact, Aj — A^ should 
be typically of order iV 1 / 6 , which can be seen heuristically as follows. 

The expected number of Aj's that are larger than 2yN — C is equal to 



E 6(A) n [2ViV-C,oo) 



arccos(w/2) , 
du 



2-c/Vn 



71 



2C 3 / 2 
37TN 1 / 4 
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for C <C y/N. Since there are exactly k parts of A that are greater or equal 
to Afc we expect that for k <C \/~N 



k ~ const ^T7i ' 

whence 

2^/N - X k ss const k 2/3 N 1/6 

which precisely means that A, — Xj should be typically of order iV 1//6 for large 
N . 

This heuristics explains the scaling in the following result (conjectured by 
Baik, Deift, and Johansson), which is the analog of Theorem^ for the edge 
u = 2 of the limit shape Q. 

Theorem 2 ([2H 03, IT7J). As N — > oo, £/ie random variables 

A, - 2V^/V 



converge in joint distribution to the Airy ensemble. 



1.9 

We now move on to plane partitions and 3D diagrams. On this set, we 
consider the following measure: 

P 3D (7r) = T^g W , qe[0,l), 
where | tt | is the volume of a 3D diagram tt and Z^jj is the normalizing factor 

^3D = ^g w = n( 1 -^ n )" n > 

it n>0 

the second equality being the classical identity due to McMahon. 

In particular, this measure gives equal weight q N to all partitions of size 
N. Is is natural to ask why we consider this essentially uniform measure and 
not some analog of the Plancherel measure on 3D diagrams (which could 
involve, for example, the number of ways to construct tt starting from the 
empty partition by adding a square at a time). The answer is that this 
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measure is the proper analog of the Plancherel measure for the following 
reason. By definition and (J2J), Ppi(A) is proportional to the number of chains 
of the form 

= \i->r> /... / AC- 1 ) / A (0) \ A (1) \ ■ • • \ \ {N) = . (11) 
such that 

A<°> = A . 

If one assembles all these A^ into an array of partitions, one sees a structure 
closely resembling a 3D diagram and, by construction, the Plancherel mea- 
sure is the distribution of the central slice of this structure induced by the 
uniform measure on all chains of the form (JTTjl . 

This also suggests that the proper way to look at 3D diagram is to separate 
them into their diagonal slices (and also suggest a common generalization of 
the measures Ppi and P 3 d to be introduced later). In particular, the diagonal 
slices of the partition from (JTJ are 

(2) -< (3, 1) -< (4, 2) -< (5, 3, 1) >- (3, 1) >- (2, 1) >- (1) , 

where A >- (i means that A and \i interlace (this relation replaces the relation 
A \ [L in (jllj) )- The proper modification of the definition (£Q) for the 3D case 
is 

©3d(7t) = {{j~ h 7Ty - ^) } C Z © ±Z , (12) 

which has the following geometric meaning. The set ©3d(tt) C M 2 is formed 
by the centers of the horizontal faces of a diagram 7r under the following 
mapping 

M 3 9 (x, y, z) ^ (y - x, z - 5±K) G R 2 . 

This is illustrated in Figure El which shows the projections of the horizontal 
faces of the 3D diagram from Figure ^ 

1.10 

The analog of the iV — > oo limit for the Plancherel measure is the q — > 1 limit 
of P3D- Indeed, as q — > 1, larger and larger diagrams become important. It 
can be shown (see e.g. jSHj) that as q — » 1 

r 3 |7r| -> 2C(3) 
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Figure 5: Projections of horizontal faces of the diagram from Figure 



in probability, where 

r = — In q . 

is a convenient parameter that goes to +0 as q — > 1. 

Thus, it is very natural to ask what is the asymptotics of the correlation 
functions 

P 3D (X C & 3D (n)) , X=(ti,hi)cZ®lZ, 
as q — > 1 and the set X is allowed to vary in such a way that 

r(U,hi) -> (u,u) G R 2 , 

while all differences 

have a limit, finite or infinite. 

At this point, it should not come as a big surprise that there is an analog 
of Theorem^ for 3D diagrams that answers this question. To state it, we 
need the following 2-dimensional analog of the sine kernel. One way to write 
the sine kernel is like this: 

K, ta(M)= ^ = i r±L. 

ixk 2m J e -i<j> w K+1 
A 2-dimensional analog of this is the following incomplete beta kernel 

K bcta (t, h-z) = ± j\l - wY -J^ , (13) 
14 



where z G C is a parameter that replaces the density parameter G M. The 
path of the integration in (J13)) crosses the segment (0, 1) for t > and the 
segment (— oo, 0) for t < 0. 

Now the analog of Theorem ^ for 3D diagrams is the following. By sym- 
metry, we can assume, without loss of generality that u > 0. 

Theorem 3 ([29j). For X as above we have 

P 3D (X C 6 3 d(vt)) -> det [K bcta (t 4 - tj, h - ^; *)1 , (14) 

J i,j=l...n 

as q — > 1, where the parameter z is the point of intersection of two circles 

{ z , z} = {\z\= e~ u/2 } f]{\z-l\ = e - u/4 ~ v/2 } , > , 

or else, in case these circles do not intersect, z = z is the point of {\z\ = 
e -«/2j w hich is the closest to {\z — 1| = e~ u / 4 ~ v / 2 }. 

1.11 

In particular, the limit of the 1-point correlation function, that is, the limiting 
density of the horizontal faces, is given by 



K be ta(0, 0; z) = = -arccos (cosh ^ - ^— J 

7T 7T \ I L J 



(15) 



Here the second equality is literally true if the two circles in Theorem El 
intersect. Otherwise, the arccosine has to be extended by or tt which, just 
like for the Plancherel measure, corresponds to trivial correlations. In other 
words, when the two circles in Theorem El do not intersect it means that we 
are trying to measure correlations outside of the limit shape. 

This limit shape, the existence of which was first established by Vershik 
[3"Tj and which was described explicitly by Cerf and Kenyon in [TUj . can be 
reconstructed by integrating the density (|T5|). This integral can be evaluated 
in terms of the dilogarithm function. The plot of the limit shape can be seen 
in Figure El 



2 Schur process and its correlations 
2.1 

It is difficult to deny that the best way to solve an asymptotic problem is to 
have a nice exact formula valid before the limit, which would make passing to 
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Figure 6: The limit shape for 3D diagrams 



the limit conceptually clear, if not completely straightforward. A wonderful 
feature of our problems is that such exact formulas are indeed available. 

In fact, exact formulas for correlation functions are available for much 
more general class of measures. These measures, which we call Schur process, 
will be described presently. By definition |2H|, the Schur process is a measure 
on sequences of partitions {A(£)}, t G Z, such that X(t) =0 for all but finitely 
many t. The probability of a configuration {A(i)} is 

PscUW)}) = 1 H S A(m± i )/A(mT l)(A m ), (16) 

where s\/^ is a skew Schur function and A m is some set of variables. The 
factor Z in (|T6|) in introduced, as usual, for normalization, that is, to make 
Pschur a probability measure. The alphabets A m and the choice of signs in 
fllfij) are the parameters of the process. In our examples, the choice of signs 
will always be understood to be such that 

■ ■ ■ C A(-2) c A(-l) C A(0) D A(l) D A(2) D ... . 
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In other words, we would like our partitions X(t) to grow up to the middle 
and then decay, just like in The general formalism of Schur processes, 
however, works equally well for any choice of signs. The above definition of 
a Schur process is slightly different from the one given in [2H], which avoids 
the awkward necessity of choosing signs, but is essentially equivalent. 



2.2 

For example, let 

A m = {q\ m \}, meZ + i, 

that is, let the alphabet A m consist of a single variable g' m ', where q is a 
parameter (e.g. a complex number such that \q\ < 1). It is a basic fact that 



s\/r(t,0,0, 



ti A ^i, Xy/x, 

, otherwise . 



where A y ii means that A and /i interlace, which can be also expressed as 
X//J, being a horizontal strip. Therefore, in this case the Schur process is 
supported on the sequences of the form 

< A(-2) -< A(-l) -< A(0) y A(l) y A(2) y ... 

and the probability of such a configuration is proportional to 

_£ m6 z+l M|A(m±i)/A(m=F|)| _ V |A(w)| 
H 2 — H ■ 

We recognize the measure P 3 d, where the diagram 7r is given by the sequence 
{A(m)} of its diagonal slices. 

The general formalism of [2H1 EH] , which will be explained on an example 
later, gives the following exact formula for the correlation functions of the 
measure P3D. We have 



P 3D (X C 6 3D (vr)) = det K 3D ((t;, hi), (tj, hj); q) , (17) 

J i,j=X...n 

where the kernel K 3 d is given by the following contour integral 



^3T>{{ti,hi), (tj,hj)',q) = 

Iff 1 $3D(ti,z) dzdw 



(2vri) 2 J\ zl=1±t J\ w \=i Tt z-w $ 3 d(*2, w) z h ^+^+\ w - h ^ 



(18) 
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In this integral, one choses the top sign if t% > t-i and the bottom sign 
otherwise. The function $3^ here is the following product 

, /, v rim>max(0,-t)(^ ~ ? I Z ) ™ , i 

^ D (t,z) = — (1 _ amz) . ^ez + i. 

At first, the integral in (|T8j) looks a little scary, but in fact it is not and, 
in particular, its asymptotic analysis, leading to Theorem El is more or less 
straightforward saddle point analysis, see [29J and below. 



2.3 

As another example, take the following specialization of the Schur process. 
Let A m be empty for m 7^ ±| and for m — ±| take the following specializa- 
tion of the algebra of the symmetric functions 

Pi^ V^, Pfc^O, k > 1 , 
where £ > is a parameter. It is well known that this has the following effect 

dim A 



|A|/2 



A ! 



Thus, in this specialization, the Schur process becomes a measure on se- 
quences 

■■•C0C0CAd0D0D... 

where the distribution of the unique nontrivial partition in the middle is 
given by 

Ppp(A) = ^£ |A| f^Y, Z PP = et. (19) 



jpp 



V |A|! 



This measure is the so-called Poissonization of the Plancherel measure, that 
is, the mixture of Plancherel measures on partitions of N for different values 
of N by means of the Poisson distribution on iV with parameter £ 

1 e N 

Ppoisson(iV) = ^ ^ , N = 0,1,2,... . (20) 

In particular, from (|20[) we can infer the value of the normalization constant 
in (HHJ). 
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It is a well known principle (known under various names such as e.g. 
equivalence of ensembles) that the N — > oo asymptotics of the Plancherel 
measures should be equivalent to the £ —>■ oo limit of their Poissonizations 
(|19|). Essentially, this is because the distribution (|20JI has mean £ and stan- 
dard deviation so for £ ~ iV — > oo it is concentrated around N. 



2.4 



The general formula of j28J |29] specializes to the following exact formula for 
the correlations of the measure Ppp (which is a limit case of the main formula 
of |E] and can be found as stated in [Sj and [T7j ) 



i,j=l...n 



Ppp (X C 6(A)) = det [KBesselte^lO 

where IK Be ssei is the discrete Bessel kernel given by 

1 ff eVZiz-z- 1 -w+w- 1 ) dzdw 



(21) 



K 



Bessel 



(2ttz) 



| Z |>H 



z — w 



1 -v+h 



Z 1 2 W 



x-y 



(22) 
(23) 



Here J n = J n (2y/£) is the Bessel function of integral order n (recall that 
x,y G Z + |) defined as the coefficient in the following Laurent series 



C2 



(z — 2' 



x z' 



The equivalence of (|2*2"Jl and (|2*3j) can be seen as follows. We have 

1 1 / 9 9 A 1 

£7; 1" W- h 1 



W -»T2 x — y \ dz dw j r T i w 
which we can substitute in (|2*2"Jl and then integrate by parts using 











oz aw 



This yields 



7 V£(z— z 1 —w+w x ) 



Z — W 



2//' 



,^(2-2 1 — w+iu 1 ) 
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2.5 

We will now explain, on the example of the formula (j21j) . the general formal- 
ism for obtaining such formulas developed in |2*H1 l2H] . Consider the following 
operator in the algebra of symmetric functions 

Since 



we conclude that 

A l A l ! 

Consider the following projection operator 



m . >a • -V: 6(A) 
, otherwise . 



Then, by construction 



Ppp(Xc6(A))= ^jf^ffi , (26) 

where we use the standard Hall inner product on the algebra of symmetric 
functions 

and where 



is the adjoint of T with respect to this inner product. The relation 

d 

dpi 
implies that 

rr = e 5 rr, (27) 

which together with 

r* • 1 = 1 (28) 
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gives another computation of the normalization constant 

Z PP = (r*r- 1,1) = . 

Our goal is to apply a similar trick to evaluate the numerator in (|2*U|) . For this 
we need a good realization of the operator ()25j) . This is naturally provided 
by the fermionic Fock space formalism (see, e.g. [20]) which will be discussed 
presently. 



2.6 



Let the vector space V be spanned by e^, k G Z + |. The infinite wedge 
space A^"V" is, by definition, spanned by vectors 

w s = e Sl A e S2 A e S3 A . . . , 

where S = {si >s 2 >...}cZ + |is such a subset that both sets 

S + = S\(Z< -§) , S_ = (Z< -i)\5 

are finite. We equip A^V with the inner product in which the basis {vs} is 
orthonormal. In particular, we have the vectors 

where A is a partition. The vector 

v% = 

is called the vacuum vector. 

By construction, the C-linear map 

s\ -> V\ 

is an isometric embedding of the algebra of symmetric functions into ArV . 
It follows from (|24j) that the operator of multiplication by pi is given in this 
realization by the natural action on A~V of the following matrix 



e x i Ae A2 a Ae A s A 



e_i A e_3 A e^s A 

2 2 2 



Pi 



a 



[ 5 M-i] fcJ 



1 
1 
1 
1 
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As we shall see in a moment, the operator (J25|) is given in this realization by 
something equally natural. 

Consider the action of the matrix unit 



En — [5fc,t5{,t] 



e sKV) 



It is clear that 



and hence 



En v x 



v x , ie&(\) 
, otherwise 



5x=\[E x 



xex 

This is already very nice, but we can actually do better by introducing the 
fermionic operators ipi which are, in some sense, the square roots of the 
operators En. 

By definition, the operator ipi is the wedging with the vector e, 

■0i v = ej A v , 

and the operator ifj* is the adjoint operator with respect to the inner product 
on A^V. One checks from definitions that 

En = ipi ip* , 

and hence the formula becomes 

1 



Ppp (X C 6(A)) = - I e^ a * H ^ r x e^ a ■ v h v d 



(29) 



2.7 



Now we want to exploit (J27)l and ()28|) . By definition, set 

^ x = Ad (V^* e -^ a ) • i/j x , 



(30) 
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and similarly define where Ad denotes the action by conjugation. Note 
that since the operators a and a* commute up to a central element, the 
ordering of these operators inside the Ad sign is immaterial. Combining 
(j27]l . ipBj) . and pty with (J2HJ) gives the following 

P PP (XC6(A))= (n**K ' V0) • ( 31 ) 

It is now a pleasant combinatorial exercise, known as the Wick lemma, to 
verify that 

II -V9,vA =det [(***J-«0,U0)] iB)W6X • 

This reduces the verification of (|21|) to proving that 

(tf x * ■ U0, «b) = K BeS sei(a;, y; . (32) 



Since fl22j) gives, essentially, the generation function for the kernel K^ssei, it 
is natural to introduce similar generation functions 



ip(z) = J2 



One then checks directly that 

= Ad (e^ a * e _ ^ Q ) ■ ip{z) = e ^ {z - z ~ r) ij){z) 
which together with 



(V>0z)V*W -v^vq) = Yl z k w k = ^^, \z\>\w\ 

£. 1 _3 _5 Z W 

2' 2' 2'--- 

implies that 

z — w 

which proves (|32j) and completes the proof of (J21j) . 
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3 Asymptotics 



3.1 

The integral representations of the type ("T""") and (J""""") are particularly suitable 
for asymptotic analysis by the classical method of steepest descent (see e.g. 
[ffl """""])■ The essence °f this method is that the asymptotics of the integral 

i MS{x) dx, M -> oo , (33) 

for a real function S(x) is dominated by arbitrarily small neighborhoods of 
the points where the function S(x) takes its maximal value. If Xq is such a 
point then S'(x ) = and hence 

S(x) = S(x ) + ^^{x - x ) 2 + 0((x - x ) 3 ) • 

Assume that S"(x ) ^ (and hence S"(x ) < 0) and introduce a new t 
variable by 

t 

X — Xq 



^-MS"{x Q ) 

Then for fixed t 



t 2 



MS(x) = MS(x ) - — + o(l) 
and we see that the leading contribution of xq to ([3*3)) is given by 



N „MS{X Q ) 



e MS(* ) I e -tV2 dx = < 



yJ-MS"( X0 ) ' 



Similar considerations apply in the case when the maximum of S(x) at x$ is 
degenerate, that is, when S"(x ) = 0, but, for example, S IV (x ) < 0. 
For similar integrals in the complex plane 



oo 



where 7 is some contour, one first deforms the contour 7 to make is pass 
through the critical points z of S(z) (also known as saddle points) in such 
a way that the real part of S(z) has a local maximum along 7 at zq. After 
that, one proceeds as with the integral (jUSJ). 
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A qualitative conclusion from this argument is that such integrals never 
have finite limit as M — > oo. In other words, whenever one expects some 
finite limit, such as, for example, we expect © to be the £ — ► oo asymptotics 
of (J22J), the contributions of saddle points computed by the method above 
will be simply zero. This makes one wonder how one can ever get some finite 
limit, because, after all, © is the £ — ► oo asymptotics of (J2*2*|) as we will see 
in a moment. The answer is that, as we deform the contours as required by 
the steepest descent method, we will be picking up residues, and those will 
have a nontrivial finite limit. 

3.2 

We now illustrate how this works on a particular example. Namely, let us 
prove that as 

xy , 
^^(-2,2), 

in such a way that x — y remains fixed, we have 

K Besse i(x,y;£) -> K sin (x - y\ 4>) , (f> = arccos(w/2) . (34) 

We have 

JK B e SS el(£, y\ = 

i rr ex P (v 7 ^ [ s g) - g ( w ; jfe)]) ^ 

JJ\z\>\w\ Z-W y/zU)' 



(2m) 2 
where 

S , (2; m) = z — z -1 — m In z . 
The critical points of the function S(z, u) are the roots of 

zS'(z; u) = z + z~ x — u (36) 

which for u G (—2, 2) are clearly given by 

z = e ±i «. 

It is easy to see that the function ?R,S(z; u) vanishes on the unit circle \z\ = 1 
and the direction of its gradient, plotted in Figure[3 is easy to infer from (|36|). 
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w 



Figure 8: Deformation of the contours for the integral (|35|) 



Thus, by the method of steepest descent we have to deform the contours of 
integration as shown in Figure|Hl after which the integral, as explained above, 
will go to as £ —>■ oo. However, as we deform the contours, we will pick 
up the residue at z = w, as is also shown in Figure |H1 It is obvious that the 

residue of the integrand in (1351) at z = w is equal to — . Thus, we 

conclude that 



Vessel 

as was to be shown. 



i f e% dw „ . 

x, V,c)— > / — = K S i n (x — y: 
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3.3 



We now turn to the asymptotics on the edge of the limit shape which for 
the Plancherel measure is described in Theorem |21 Recall that near the edge 
u = 2 of the limit shape the distance between two consecutive "downs" is 
of order iV 1 / 6 and hence also the probability to observe a "down" in any 
particular position is of order iV -1 / 6 , where N is the size of the partition A. 
For the Poissonized Plancherel measure Ppp the expected size of A is £. From 
Theorem |21 we, therefore, expect that 

e /6 KBessei (2^ + x^ 6 , 2^ + y^ 6 ) - K Airy (x, y), £ - oo . (37) 

We will now use the steepest descent method to deduce (J37)l from The 
limit ()37|) is the key ingredient of the proofs of Theorem El given in 0GZ]. The 
difference between (|37|) and Theorem El is a technical, but rather standard, 
collection of estimates that allows one to deduce the convergence of joint 
distributions for Plancherel measures from the convergence of correlation 
functions of their Poissonizations. 

3.4 

We now start the asymptotic analysis of the integral (f33|) in the u = 2 regime. 
What happens in this case is that the two nondegenerate critical points e ±l< ^ 
of S(z, u) coalesce at z — 1 to form one degenerate critical point 

S{z,2) = \{z-lf + 0{{z-lY). (38) 

Since S(z,u) appears in the exponent in (f33j) . this starts looking like the 
Airy integral pOjl . and, indeed, this is precisely the source of the appearance 
of the Airy function in ([37)1 . In fact, this is the fundamental reason why the 
Airy function is as ubiquitous in asymptotics as the sine function. The sine 
function appearance in the steepest descent asymptotics is usually produced 
by two complex conjugate critical points. When the parameters are tuned 
so that these two points coalesce, the Airy function enters the scene. This is 
all very standard and is at length discussed, for example, in [1U] and many 
other places. 

The integral (fTUjl converges, but only conditionally, due to ever faster 
oscillation of the cosine function for large s. A much better integral can be 
obtained from by shifting the contour of integration as shown in Figure 
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Figure 9: A better contour for the Airy integral (jlUJI 



El The contour in Figure 01 goes along 2 of the 3 downslopes of the graph 
of 3ft(zs 3 ), plotted in Figure^] and commonly known as the monkey saddle. 
The function e ts3 ^ is very rapidly decaying on this contour. 




Figure 10: Plot of K(is 3 ) 



3.5 

It is clear that for u = 2 all the action in the integral ()35|) is happening in 
the neighborhood of the critical point z — w — 1, so we want to zoom in on 
this point. We also want the integration contour in z to pass through the 
downslopes of the graph of the real part of (J38J) . while the w-contour should 
go along the ridges of the same graph (because the of the opposite sign of 

w 'i ~7? ) ^ erm i n Q35j)), We thus deform the integration as shown in the 
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left half of Figure ^2 and introduce the following new integration variables 



. z' , vJ 

Z = 1 + 1—TJ7 , W = 1 + I 



£ 1/6 



£1/6" 



(39) 



For large £, the integration contour in z' and w' will look like the contour in 
the right half of Figure The scaling in ()39j) is chosen so that to make the 




7 IV 




IZ 



Figure 11: Contours for the edge (u = 2) asymptotics 
exponent in 

^S(z,2) = - l -(z'f + o(l) 

stay finite as £ — > oo. 

More generally, the leading order expansion of S(z, u) near (z, u) 
with respect to both arguments is given by 

(z,2 + ^ J = -ix*' - % - (z'f + 0(1) , £ ^ oo . 

Plugging this into the integral (|35|) we obtain the following 



e 1/6 Kbcssci (2 + < 1/6 , 2 + y£ 1/6 



(27T) 



e -i(^') 3 /3-«z'^+j(«'') 3 /3+ra'l/ 



tteW, (40) 



where the integration contour (rotated by 90°) is shown in the right half of 
Figure HU 
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3.6 

It remains to show that 

1 f f e - ^ 2 '-* 3 / s ~ iz ' x + i ( w ') 3 /3+iw'y 



(2?r) 2 J J %{z' - w') 

Ai(x) Ai'(y) -M{x) Ai(y) 



dz'dw' = 

(41) 



x-y 

This is entirely parallel to the argument in Section 12.41 (see the more general 
discussion in Section 2.2.4 of [2H])- We have 

i (3 (9 \ 

^—iz'x+iw'y I | j ^—iz'x+iw'y 



x — y \ dz' dw' t 
Inserting this into the LHS of (|41j) and integrating by parts we obtain 

e -i(z') 3 /3-iz'x+i(w') [i /3+iw'y 

dz'dw' = 



(2tt) 2 J J i(z' - w') 

— , [[(Z' + W') e -i^-iz'*+iH£+™'v dz l dw > . 

(2 7 r) 2 (x- Z /)yy 1 ; 

The formula (fTUj) translates into 

— / e~ l ^^ 3 ~ lz ' x dz' = Ai(x), 
2n J 

because our contour of integration in z' (see Figure ITT|l is the negative of the 
contour in Figure 01 It is also clear that 

— / w'eW Mv dw' = Ai'(y). 
2n J 

Thus, equation (|4*Tj) follows. 



3.7 

After finishing these computations, it may be useful to look at them from a 
more abstract point of view. Above, we tried to be absolutely concrete and 
explicit, but in fact we used only very general properties of the correlation 
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kernel (|35|). Namely, all that was needed was the information about the lo- 
cation of the critical points of the function S in (|35j) and the residue of the 
integrand in ()35|) at z = w. The fact that the two critical were complex con- 
jugate, which was important for obtaining the sine function, is an automatic 
consequence of the fact the integrand is real for real values of the argument. 
Even less was required for the Airy function asymptotics, the only essential 
ingredient in which was the coalescence of the two critical points. 

This means the same method is applicable to much a wider variety of 
problems producing the same or equivalent results. This phenomenon, when 
the final asymptotic result proves to be insensitive to the fine details of the 
original problem, is known as universality. The simplest example of this 
phenomenon may be the central limit theorem which says that, under very 
general assumptions, sums of independent, identically distributed, but other- 
wise arbitrary random variables converge, scaled and centered, to the normal 
distribution. It is widely believed that the universality principle applies very 
generally. From a mathematical perspective it is often quite mysterious why 
this should be the case, but sometimes, like for the central limit theorem, 
there is a transparent explanation. Similarly, our exact formulas for the cor- 
relation functions of the Schur process, combined with the steepest descent 
argument, offer a simple explanation for the appearance of the sine and Airy 
kernel. 
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